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Abstract 

The depth rule is a level truncation of tensor product coefficients expected to be sufficient 
for the evaluation of fusion coefficients. We reformulate the depth rule in a precise way, 
and show how, in principle, it can be used to calculate fusion coefficients. However, we 
argue that the computation of the depth itself, in terms of which the constraints on tensor 
product coefficients is formulated, is problematic. Indeed, the elements of the basis of 
states convenient for calculating tensor product coefficients do not have a well-defined 
depth! We proceed by showing how one can calculate the depth in an 'approximate' way 
and derive accurate lower bounds for the minimum level at which a coupling appears. It 
turns out that this method yields exact results for su(3) and constitutes an efficient and 
simple algorithm for computing su(3) fusion coefficients. 
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1. Introduction 

By now various methods have been proposed to calculate fusion rules in WZNW 
models [1-7]. The first method to be proposed was the depth rule of Gepner and Witten 
[1] . However it has not received much attention since it did not appear to be very practical, 
mainly because of the difficulties involved in calculating the depth. Therefore the depth 
rule has not been tested seriously as a calculation tool for fusion coefficients. Indeed, 
except for su(2) : only low level examples of su(3) [1] and E 8 [8] have been worked out. 
In this work we examine carefully the depth rule as a computational tool for fusion rules. 
Section 2 contains a review of the original derivation of the rule (see also [9]). 

Recall that a fusion rule tells us which primary fields can arise in the OPE of two 
given primary fields. More abstractly it defines a product x as 

where A, \i and v are primary fields. The fusion coefficient N^J v gives the number of 
times v appears in the product A x \i. Notice further that N^ k J v = N^ k J Cu , where Cv is 
the field conjugate to v. For WZNW models the spectrum generating algebra is a Kac- 
Moody algebra g at level k and the primary fields are in one-to-one correspondence with 
integrable representations of g [1,10]. 



2. The depth rule of Gepner and Witten. 

Let us start by writing explicitly the commutation relations for the untwisted Kac- 
Moody algebra g at level k [11]: 

[H l m , Hi] = km5 lJ 5 n+m , 



a 

n-\-m 



[Hi n ,E%]=a i E : 

[E^Ei] = ^(a-H n+rn + kmd n+m , ) if a = -(3 (2.2) 



= e(a, (3) if a + (3 is a root 

= otherwise 

Here a and f3 denote roots and e(ct, 0) is a cocycle. To each node (i.e. each simple root) 
of the extended Dynkin diagram one can associate at least two distinct sit (2) subalgebras. 
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In particular for the zeroth root, whose finite part is —9 (9 being the longest root), the 
explicit forms of the generators and J3 of the two su(2) subalgebras we will use are 

(i) Eq , Eq B , 9 ■ H Q = [Eq,Eq 9 ] 

(ii) E~ e , E e _ x , k-9-H = [E- 9 ,E d _ l ] 2 ' 3 

The finite Lie algebra g associated with g is recovered from (2.2) by setting m = n = 
(this corresponds to the homogeneous gradation). Hence the first su(2) subalgebra in 
(2.3) is generated by finite Lie algebra operators while the second one has an intrinsic 
affine structure. Now let A denote an affine weight of positive level k, whose finite part A 
is an integrable highest weight. There is a one-to-one correspondence between such affine 
weights and WZNW primary fields, and we denote them by the same symbol. To every 
field A we associate a state A(0)|0) = |A). The degenerate nature of the WZNW primary 
fields translates into the following constraints for the singular vectors [1]: 

(£-*)(*.*)+i|A) = (2.4) 

(^i) fc-(A ' fl)+1 |A) = (2.5) 

Of course there are also singular vectors with 9 replaced by an arbitrary simple root a. 

These singular vectors lead to constraints on three-point functions. The constraints 
arising from the first type of singular vectors (2.4) being associated with finite Lie algebra 
generators, they are taken into account in ordinary tensor product coefficients of finite 
Lie algebras [12]. On the other hand, the other singular vectors yield new constraints 
on the ordinary tensor product coefficients which depend explicitly on the level k. Since 
(A, 9) > (A, a) for any root a of g, the maximal constraint is obtained by considering 
the longest root, i.e., precisely from (2.5). Let us then derive the explicit form of this 
constraint. For this consider the three-point function 

(((E e _ i r\)(z)^(z i y'(z 2 )) (2.6) 

where \i" is associated to a weight (also at level k) whose finite part jl" belongs to the 
representation with highest weight jl (similarly for u"). This vanishes if p > k — (A, 9) as 
a result of (2.5). By considering the E B n {z) as the modes of the field E d (Q in a Laurent 
expansion about z, i.e. 

E e (o = J2 E n(z)((-zr n - 1 

n 
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one can write 




(Factors of (2ni) 1 in front of contour integrals are everywhere implicit). The three-point 
function (2.6) can thus be transformed into 

Let us now deform the integration contour of every variable Q such that it circles around 
the singular points z\ and z 2 . Expanding the fields E 9 (Q) in terms of their modes and 
performing the integrations (in which only the mode m = contributes) one finally gets 

t t) , {z ^ {z _ ^ w»)[(Wfe)[(^M> 

As already mentioned, this vanishes if p > k — (A, 9). Since all the terms in the sum are 
independent, they must all vanish separately, i.e. 

(A(z)[(Sj)y / ](z 1 )[(Sj)"-'i/ / ](z 2 )> = if p>k-{\9) (2.7) 

This holds for any value < I < p. 
We conclude that if 

(Eft V = /*V o 

(2.8) 

for li + l 2 = p > k — (A, 9), the nontrivial three-point function [Xji'v') must vanish. Here 
fi' (reps, is') is characterized by the fact that one can act l\ (reps. 1 2) times on it with 
Eq 9 without leaving the representation, i.e., one obtains a non-zero field fi" (resp. v"). 

Let us now define the depth of the field ji' as [1] 

= max(I) such that (E~ e ) l /j' ^ (2.9) 

and similarly for v' (we will come back to this definition later and make it more precise). 
Since l\ + I2 < gL' + d u i we conclude that 



(X(z)ii'(z i y(z2)) = if k <d», + (\,9) 

3 



(2.10) 



This is the explicit form of the constraints we were looking for. Of course, irrespective of 
the value of k, it is always true that 



(X(z)^'(z i y(z 2 ))=0 if N- Xpp = 

where denotes the coefficient associated with the finite tensor product A <8> ft <8> V D 0. 

Since 

= N- X f 

where Cv is the conjugate of 9, N Xjlj? is equal to the number of times Cv appears in the 
tensor product A <g) ft. (Recall that a necessary condition for Nxfa? 7^ is the existence of 
a ft' (resp. v') in the highest weight representation ft (resp. v) such that A + ft' + v' = 0). 
We summarize these results in the following depth rule [1]: 

The coupling (A/iV) vanishes if either k < +d L ,> + (A, 9) or N\fii? = 
0. 

Before proceeding to a detailed analysis of the depth rule, let us emphasize the fol- 
lowing remarkable fact. As is clear from (2.9), the depth d^> is a property of the finite 
part of the weight ji' (i.e. it involves a generator of the finite Lie algebra). Hence, starting 
from a singular vector involving the action of a generator not in the finite Lie algebra, 
namely E e _ 1: which amounts to a constraint on the affine weight, one ends up with a con- 
dition whose only residual affine characteristic is the level. Thus in applying this rule, all 
computations are done with finite weights. 

Notice further that the constraint (2.5) is void when k — > 00. In that case we are left 
with purely finite Lie algebra constraints, so that 



^ N x% = N-x^ (2.11) 



3. Deepening the depth rule 

3.1. States vs. weights. 

The derivation of the depth rule, reviewed in section 2, depends crucially upon con- 
sidering correlation functions, which themselves involve states. We have pointed out that 
WZNW primary fields are in one-to-one correspondence with affine weights A of positive 
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level k whose finite parts are highest weights A. These in turn are in one-to-one corre- 
spondence with states |A). This one-to-one correspondence between weights and states no 
longer holds if A' has a finite part A' which is not a highest weight, but rather a degenerate 
weight with multiplicity mult (A') > 1. To this weight correspond mult (A') states |A^), 
i = l,..., mult (A'). This clarification is crucial since the depth is a property of states, not 
of weights, and it is the states which characterize a given coupling. The expression (2.9) 
for the depth should then be sharpened: 



= max(l) such that (£ "")V t) ) / (3.2) 

Similarly the notation (A/iV) should now stand for the correlation (Xfi'^v'^) where the 
doublet (j, £) can take N\fii> values. 

An aesthetic drawback of the original formulation of the depth rule is the lack of 
symmetry between the three states. This can be remedied by considering, instead of a 
three-point function, the following four-point function 

which vanishes when p > k + 1. Here / denotes the identity field, associated with a null 
finite weight. Following the steps of the previous section, we obtain that the three-point 
function (A' (i) (z 1 )^(z 2 )i''^(z 3 )) vanishes if: 

k < d u > + d v i + dy (3.3) 

where p! + v' + X' = 0. Notice that this reduces to (2.10) when A' = A, since d\ = (A, 6>)). 
Eq. (3.3) is a symmetrized version of the depth constraint, which makes manifest the 
symmetry of the fusion coefficient N^ x under the interchange of any two labels. 

However, for practical calculations, it is preferable to fix two of the states, that is to 
set A^ = A and v'^ = —Cv and consider then the three-point function 

(X(z 1 )^ l) (z 2 )(-Cu)(z 3 )}, i = l,...,iV^ (3.4) 

Here we stress that the /tt'^'s are understood to be the states which, when coupled to A 
and —CP, give the scalar representation. The advantage of (3.4) is that only one out of 



the three weights is neither a highest nor a lowest weight, which greatly simplifies the 
calculations. Since the depth of —Cv is obviously zero, (3.4) vanishes whenever 



k<d^. +(A,0) . (3.5) 



5.^. Threshold level k^ and the reformulation of the depth rule 

When k — > oo, there are no affine constraints, so that (\fi'^(—Cv)) ^ 0. This actually 
holds true for all values of k sufficiently large that one does not 'hit' the singular vector 
(2.5). For k < d^ + (A, 6), the coupling must vanish as a result of the depth constraint. 
As a result, the smallest value of the level at which the coupling will be nonzero is 

= cty o + (A, 9) . (3.6) 

For each of the N\ni> couplings associated to the tensor product A <8> ft <8> V D 0, there 
is then a threshold level, Kq \ such that for /c > Kq the coupling is non zero, while for 
k < &q it vanishes. Accordingly, the fusion coefficients are completely determined by the 
tensor product coefficients and by the possible values of fcp given above. Suppose now 
that the values of k^ are ordered such that k^ < kQ +1 \ We can then reformulate the 
depth rule in a precise way as: 

N$® v = max(j) such that k > k^ ] and N~ XflP ^ 

= if k < k ( Q } or N-^ = ^ 

with k^ defined in (3.6) in terms of the depth. 

It remains to see how to calculate the depth. But first we introduce an auxiliary 
concept, the depth charge. 



3.3. Introducing the "depth charge" ag. 

Suppose we have a basis of states such that applied to an element of this basis 
either vanishes or gives another single element. That is, one never gets linear combinations 
of basis elements by applying E^ 6 '. Then the depth is essentially the maximum number 
of times one can act on a given state with E^ 6 without leaving the representation. If 
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we denote by \jm) the projection of this state on the su(2) subalgebra associated to 9 
( J 2 = j(j • + 1) and m is the J 3 eigenvalue) then the depth is simply j + m. We will write 

a = 2j (3.8) 

and call ag the depth charge. For a state \^'^), m is the eigenvalue of 9 ■ H , and is given 
by 

(fi',9) = 2m 
Therefore, the depth can be reexpressed as 

2cV w = a,(p' (0 ) + (/*',*) (3.9) 

The calculation of then boils down to that of ag. In terms of ag, (3.6) reads 

2A#> = ae (»{ i) ) + (\ + v,9) (3.10) 

where we used (/t/, 9) = {Cv — A, 9) = {v — A, 9). 

3.4- Calculation of the depth. 

In order to test the depth rule, one needs to know how to calculate the depth charge. 
For this one must answer the following three questions: 

1. How can we describe states explicitly? 

2. How can we characterize the three states appearing in a three-point function? 

3. Given a state, how do we calculate the depth charge ag! 

These questions are looked at in the next section. From now on the discussion will be 
restricted to su(N) for simplicity. However the results presented here can be generalized 
to other algebras. 

4. The depth machinery 

4-1. A basis for states: standard tableaux. 

An su(N) integrable highest weight A = J^iLi 1 = C^i> • • • > -V/v-i), where oo l is 
the ith fundamental weight, has a natural representation in terms of a Young tableau with 
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Ai + . . . + Atv_i boxes on the first row, A2 + . . . + Xn-i boxes on the second row, up to 
Aat_i boxes on the (N — l)th row. Standard tableaux (also called semi-standard tableaux 
[13]) are Young tableaux in which each box is numbered according to the following rule. 
Let Cij be the number appearing in the box on the ith row and the jth column. These 
numbers must satisfy the following constraints: 



l<Cij<N , 



< Ci+1, 



In words, the numbers are non-decreasing from left to right and strictly increasing from 
top to bottom [14]. 

The allowed standard tableaux of shape A are in one-to-one correspondence with the 
states in the highest weight representation A. Weights A' with multiplicity mult(A') are 
associated with mult (A') different standard tableaux. Thus standard tableaux provide a 
convenient basis for states. 

Notice that the weight of a standard tableau is obtained by adding the weights of 
all its boxes. The latter are easily obtained from the fundamental representations: u k is 
associated with a single column of k boxes containing the numbers 1,2, ... ,k. 

Example 1: 

For su(3) the standard tableaux of the highest weight states of the fundamental rep- 
resentations are 



(1,0) ~ a , (0,1) 



(0,0) 



from which we conclude that (—1, 1) <-» [2] and (0, —1) <-> The adjoint representation 
(1,1) contains 8 standard tableaux: 



1 


1 


1 


2 


1 


3 




1 


1 


1 


2 


1 


3^ 


2 


2 


2 


3 


2 




2 




2 






_3 




_3 




3^ 




3^ 




3 




;m) (■ 


-1,2) (0,0) 


(2, 


-1) (0,0) (1, 


-2) (■ 


-2,1) (- 


1, 


-1) 



Notice that to the doubly degenerate weight (0, 0) one associates two distinct standard 
tableaux, i.e. two distinct states. 

In the following we will use an equivalent representation of this basis, the so-called 
Gelfand-Tsetlin patterns [15]. To a given standard tableau we associate the following 



triangular array of numbers: 



Pi 



(N-l) 



3 {N) 
Pn 



(N-1) 



p { ? p { ? 



(4.2) 



Pi 



(i) 



such that (iff' is the number of boxes containing numbers less or equal to j in the ith row 
of the standard tableau . 

Example 2: 

The following standard tableau and Gelfand-Tsetlin (GT) pattern corresponding to 
the weight (—2, 1, 0) in the representation (1, 2, 1) of sit(4) are equivalent: 



4 3 10 
3 2 1 
3 1 
1 



1 


2 


2 


4 


2 


3 


4 




3 









4-2. Identifying states in tensor products. 

Given a triple product involving a highest weight state A and a lowest weight state 
—CP, what are the standard tableaux of shape Jjl and weight Jjl' = Cu — X that can appear in 
the product? To answer this question we first write down the mult(/i / ) standard tableaux of 
weight Jjl' . Among these, the standard tableaux contributing to the product A®//(g>z/ 
are those which, when inserted column by column from right to left into the Young tableau 
of A, still give a regular tableau at every step. The insertion is made by adding to the ith 
row of the Young tableau the boxes of the standard tableau marked by i [16]. 

Example 3: 

Consider the product (1,1)<8>(2,1)d(1,0). We are looking for standard tableaux of 
shape (2, 1) and weight (0, —1). There are two of them 



1 


3 


3 


and 


1 


2 


3 


2_ 






_3 





Adding the first standard tableau to 



(the Young tableau of shape (1,1)) column by 



column from right to left, one ends with a non-regular tableau after 2 steps: the third 
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row has two boxes and the second row only one. Thus the first standard tableau does not 
contribute to the product, while the second one does (as is easily checked). 

We now know how to characterize the states in triple products involving one highest 
and one lowest weight state. In general it is also possible to describe all the states in a 
given triple product, but this requires the introduction of further concepts which are not 
central in our argument, and consequently they are introduced in the appendix only. 

4-3. Calculation of the depth charge for the standard tableaux. 

A GT pattern indicates the chain of embeddings u(2) C w(3) C . . . C u(N) where the 
su(2) C u{2) subalgebra is defined with respect to the root ct\, the su(3) C u(3) subalgebra 
is defined with respect to the roots a± and «2, etc [15]. As a result the spin j of a given 

state with respect to the su(2) subalgebra associated to ol\ can be read off directly from 

(2) (2) 

the GT pattern. Denoting by a ai twice this spin, one has a ai = (3\ — fj ^ , where the /3's 
are defined just after Eq.(4.2). Thus the properties of the states with respect to the su{2) 
subalgebra associated with a.\ are easily extracted. 

On the other hand the action of the operators E^ e on standard tableaux or GT 
patterns is not so simple. Acting with E^ e on a standard tableau generically yields a 
linear combination of standard tableaux with the same weight. Hence a GT pattern does 
not have a well-defined depth charge. A first consequence of this observation is that only 
lower and upper bounds for ag can be calculated for states whose weights have multiplicity 
> 1. These bounds are the minimal and maximal values that ag can take for all the 
states corresponding to a given weight (not just for the states contributing to the tensor 
product). As far as ko is concerned only the lower bound is relevant, but does not turn 
out to be very useful. 

Thus one faces a problem with the depth rule approach for calculating fusion rules. In 
this scheme fusion coefficients are truncated tensor product coefficients. However the basis 
convenient for the calculation of tensor product coefficients (namely, standard tableaux) 
is not convenient when time comes to evaluate the degree of truncation appropriate to a 
certain level k in fusion rules. 

3 Notice however that for a coupling (Xfi'v') such that v' = —Cv and mult(/i') = N\nv we can 
say better. Indeed if we know all possible values of ae which corresponds to the weights, these 
values can be used in eq.(3.10) to calculate the possible values of ko- 
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4-4- Introducing the crystal depth charge a c B . 

In spite of this situation, we will try to make further progress, motivated by the 
following observation. It is well known that fusion rules for su(N) at level k are equivalent 
to restricted tensor products in the quantum version of su(N) with q given by [17]: 



On the other hand, at q = the structure of the quantum group simplifies considerably 
[18] (the group 'crystalizes'). For our purpose it is remarkable that an explicit realization 
of a crystal base is given by the standard tableaux [19]. Furthermore, it is a particularity 
of q — that standard tableaux have a well-defined value of ag. This prompts us to use 
these values of ag computed for the quantum group at q = 0, which we denote a B (the c 
stands for crystal), in the formula for kg. Although the logical motivation for making such 
a naive extrapolation from q = to q = exp{2ni / (k + N)} may look weak, it turns out 
that this procedure leads us to rather accurate results and, for su(3), exact results! 

One complication did arise, however. There is a one-to-one correspondence between 
the three couplings 



(This correspondence can be seen for instance via the Berenstein-Zelevinsky triangles, to 
be introduced in section 6.) We found that the values of k^ calculated in this way differed 
in general, for the 3 couplings above. It was therefore necessary to take the maximum of 
the values so obtained. We symbolise this by rewriting (3.6) as 



As already pointed out, standard tableaux or GT patterns naturally single out the 
su(2) subalgebra associated with the root a±. Thus, in order to calculate ag we will make 
a suitable 'rotation' of the GT pattern such that ag will be the difference of the new entries 
sitting on the line second from the bottom. But since this operation is justified at q = 
only, it is a c e that we obtain in this way. Unfortunately, in order to define the explicit 
operation of 'rotation' some more technology must be introduced. 




(V«)(-Ci/)>, <^)(-CA)>, (v\'m(-Cvl)) 




maxpermlc^y + (A, 6)} . 



(4.3) 
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Let us define the operators t j which act on the GT patterns as follows [20] : 



(fc) 



(0 



(3f ] if k^i 



>(0 



If Pj-i^ is absent it is set equal to oo and similarly if /?■ ^ is absent it is set equal to 
zero. Then let us introduce the operator 

a = t\ti . . . £jv_i 

for su(N). These operators have been introduced to describe the action of the Weyl group 
directly on GT patterns. The basic Weyl reflections Sj (with respect to ctj) are simply [20] 



This is equivalent to the action defined by Lascoux-Schiitzenberger directly on the standard 
tableaux [21]. 

Example 4: 

Here are various operations on a particular state of weight (0, 1, —1, 0) in the highest 
weight representation (0, 1, 1, 1) of su(5) producing in the end a state of weight (0, 0, 1, —1): 



1 


1 


3 


2 


2 


5 


4 


4 




•5 







33210 33210 



3 2 2 
3 2 
2 2 
2 



t 4 



3 3 10 
3 2 
2 2 
2 



3 3 2 1 
3 3 10 
3 2 1 
2 2 
2 



*2 



33210 33210 



3 3 10 
3 2 1 
3 1 
2 



ti 



3 3 10 
3 2 1 
3 1 
2 





1 


1 


2 




2 


3 


4 




3 


5 






•5 







This sequence of operations is nothing but the action of o. 

Now, o turns out to be the desired 'rotation' operator. Indeed, from the results of 
[20] one can show that 



_— 1 771 a i „ 

a hiQ a 



12 



Example 5: 

The action of a = t\t2 on the adjoint representation of su{3) is shown on Fig.l. One 
sees that this action amounts to a rotation of the peripheral weights by 120° and to an 
interchange of the two inner states. The 9 direction has been rotated to the ct\ direction 
except for a 'sign', which illustrates the fact that Eq (a J_ operator) is actually related 
to Eq 1 (a J + operator). 

In terms of GT patterns this means that a c e can be calculated as follows. Let us write 

GT = a(GT) 
where the entries of GT are denoted j3) ' . Then 



Example 6: 

The following state corresponds to the weight (0, —1,0) in the sit (4) highest weight 
representation (1,1,1): 



1 


3 


4 


2 


4 




3 







3210 3210 

2 11 a 3 10 
11 ~" 3 
1 2 



a% = 3 - = 3 



Similarly the value of a c d for the standard tableau considered in example 4 is 2. 



4-5. A simple rule to calculate a c B for su(3). 

In this section we derive a simple rule for computing the value of a@ directly from 
standard tableaux in the case of su(3). In this case a = t\t<ii but since t\ affects only the 
last row of the GT pattern, it is sufficient to consider only t 2 - Hence one has 



a b c a b c 
d e d e 
f f 



where 

d = a + max(6, f)—d 
e = min(6, /) + c — e 
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Dropping from the outset columns of three boxes from the standard tableau, one can set 
c = (this amounts to setting the irrelevant u(l) charge of u(l) <g> su(N) C u(N) to zero). 
Thus 

ag = a — d + e + max(6, /) — min(6, /) 
= a-d + e+ \b-f\ 
This is easily checked to be equivalent to 

a c B = ^ (number of l's — number of 3'sj (4.4) 

columns 

In other words, the value of a% of a given standard tableau can be calculated by adding 
the values of a c Q corresponding to each of its columns. Furthermore, in each column the 
value of a c e is just the number of l's minus the number of 3's (in absolute value). 

Example 7: 



a* 



1 


1 


2 


2 


2 


3 


2 


3 


3 





= 3 



5. Testing the depth rule with a c e . 

We are now in a position to test the depth rule with crystal depth charge through 
various examples. 



5.1. The su(3) case. 
Example 8: 

Let us first consider a classic example, namely the triple product (l,l)(g)(l,l)<g)(l,l). 
It is well-known that (l,l)<g>(l,l) = 2(l,l)© — In fusion rules one copy of (1,1) arises 
at level 2 while the other appears at level 3. Indicating the value of ko by a subscript one 
has 

(1, 1) <8) (1, 1) = (1, 1) 2 © (1, 1) 3 © - - - 
To the weight (0, 0) correspond two standard tableaux of shape (1, 1): 

and 
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Their values of a c e are respectively and 2 (c.f. (4.4)). Thus here (3.10) gives directly 
ko = 2, 3, which is the correct result (permutations are not necessary since the three 
weights are equal). 

Using a computer program we have tested the depth rule through (4.3) with the depth 
given by (3.9) and ag calculated at q = by (4.4). We computed all fusion coefficients up 
to level 10 and found perfect agreement with the results obtained from the Kac- Walton 
formula [3] . This led us to suspect that the depth rule with crystal depth charge is correct 
for su(3). Actually, this can be proved directly by relating our algorithm to other known 
algorithms for su(3). We will return to this in the next section. 

5.2. An sit (4) counterexample. 

Here we present an example which shows that the depth rule with crystal depth charge 
is not correct in general. The simplest counterexample that we found is the following sit (4) 
product: 

Example 9: 



(l,2,l)(g)(l,2,l) = 4(l,2,l) 5 e(l,2,l) 6 



(5.2) 



To the weight (0,0,0) in the representation (1,2,1) correspond seven standard tableaux, 
among which five (obtained by the method of section 4.2) contribute to the product: 



1 


1 


2 


3 


2 


3 


4 




4 









1 


1 


2 


3 


2 


4 


4 




3 









1 


1 


3 


3 


2 


2 


4 




4 









1 


1 


3 


4 


2 


2 


4 




3 









The corresponding values of ag, calculated according to the rule given in section 3.3, are 
respectively 0,2,2,2,4. Hence formula (3.10) yields the following values of ko'. 



ko = 4 , 5 , 5 , 5 , 6 



which disagrees with 5.2. 
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6. su(3): relation with other methods. 



In this section we show that, in the case of su(3), the depth rule with crystal depth 
charge is equivalent to all other known algorithms for computing su(3) fusion rules. This 
equivalence is obtained via the Berenstein-Zelevinsky (BZ) triangles, introduced in section 
6.1. We then reexpress ko as given by eq.(3.10) in terms of the entries of these triangles. 
Next we recall that the generating function for su(3) fusion rules obtained in [5] (and 
whose correctness was proven by Cummins [22]) can be translated into a simple formula 
for ko. This in turn can also be reexpressed in terms of the parameters of the BZ triangles. 
The resulting formula for k is equivalent to that obtained in [23] by other methods. We 
then show that it is also equivalent to the formula derived from the depth rule with a c e . 

6.1. Berenstein-Zelevinsky triangles for tensor product coefficients. 

Consider the set of three su{3) highest weights (Ai,A2), (^1,^2) and (^1,^2) (we 
give here the Dynkin labels). Berenstein and Zelevinsky [24] showed that the number of 
triangles one can construct according to the following rules: 

ai + a 2 = Ai 
CI4 CI3 + CI4 = A2 

as as such that ^ + a 5 = m 

a 2 a 6 o 6 + a 7 = \i 2 

a\ ag as a-j aj + as = v\ 

ag + a\ = v 2 

gives the value of N\^ D . Such triangles make manifest most of the symmetries of the tensor 
product coefficients. 

Example 10: 

Corresponding to the coupling (2, 2) ® (2, 2) ® (2, 2), three BZ triangles can be con- 
structed: 

1 2 

2 2 11 

2 2 ' 11 ' 

0220 1111 2002 

and accordingly the multiplicity of the coupling is 3. 
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«2 + ^3 = <26 + as 

as + 05 = ag + as (6.2) 

a 5 + a 6 = a 2 + a 9 



The states involved in a specific coupling can be read off a triangle as follows. Consider 
the product A <g> p, = Cv © . . . associated with the BZ triangle (6.2). The state \p,'^ > (of 
weight Cv — A) in this coupling is described by the GT pattern 



A*i + A*2 - a 4 \i 2 - a 6 
A*i + A*2 - a 4 - a 2 



(6.3) 



Example 11: 

The GT patterns and corresponding standard tableaux (in the representation jl) as- 
sociated with the three BZ triangles of the last example are (in the same order) 



4 2 
4 
2 



1 


1 


2 


2 


3 


3 





4 2 
3 1 
2 



1 


1 


2 


3 


2 


3 





4 2 
2 2 
2 



1 


1 


3 


3 


2 


2 





£.2. ko from BZ triangles. 

Consider the BZ triangle (6.2) describing the coupling between A, v' = —Cv and \x' 
given by (6.3). The value of ag(n') can be obtained by acting on the GT pattern (6.3) 
with a = tit 2 (actually t 2 alone is enough in this case) and by taking the difference of the 
entries on the second line. This gives 

a e(A0 = a 4 + «7 + |«5 — «2| 

Therefore one has 

agin') + {v + A, 9) = 2 max(a 4 + v\ + z/ 2 , a-? + Ai + A 2 ) 

The permuted versions of the above are evaluated similarly. As a result, one can rewrite 
(3.10) in the form (There are only 3 permutations here because interchanging v and A does 
not change the value of the crystal depth charge a^/u').) For a 2 < as < as, this boils down 
to 

k = max{a 7 + Ai + A 2 , a 4 + z/i + v 2 } ■ 
But (a 4 + ui + z/ 2 ) — (a7 + Ai + A 2 ) = as — a 2 > 0, by assumption, so we obtain 

k = a 4 + v\ + v 2 ■ (6.4) 
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It can similarly be shown that if we assume instead a<i < a§ < as, the formula (6.4) remains 
valid, and so holds for a 2 < min(as, as). 

Example 12: 

From the above formula one readily computes the values of k for the three triangles 
of example 10 to be respectively 4, 5 and 6. (Notice that for these 3 triangles the inequality 
^2 < min(a5,ag) is satisfied. For less obliging triangles, however, one simply uses their 
symmetry to rotate them until the inequality is obeyed.) 



6.3. ko from the decomposition of BZ triangles into elementary couplings. 

Recall that every coupling A <S> ft <S> v can be decomposed into a product of elementary 
couplings Ei [25]. For su(S) the 8 elementary couplings are 



E x = (1, 0)(0, 1)(0, 0) E 2 = (1, 0)(0, 0)(0, 1) E 3 = (0, 0)(1, 0)(0, 1) 



1 

1 1 

0000 1000 0100 



E 4 = (0,1)(1,0)(0,0) E 5 = (0,1)(0,0)(1,0) E 6 = (0, 0)(0, 1)(1, 0) 

1 

1 



0000 0010 0001 



E 7 = (1, 0)(1, 0)(1, 0) E 8 = (0, 1)(0, 1)(0, 1) 


1 10 
10 1 

1 1 

However there is some redundancy in such a decomposition since different products 
of elementary couplings can be equivalent. Since the BZ triangle of a product coupling is 
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simply the sum of the BZ triangles of the factors, one sees that EiE s E 5 = E 7 E 8 , since 
they both have the same BZ triangle: 



E1E3E5 = EfEg 





1 1 
1 1 
110 



For su(3) this is the only redundancy [25], and to obtain a unique decomposition for 
a general coupling one simply has to forbid one of the products E 7 Es or E\EzE§ to 
appear. This leads directly to the construction of a generating function for tensor product 
coefficients. 

It was conjectured in [5] that there is at least one choice of forbidden couplings which 
allows for the construction in a simple way of a generating function for fusion coefficients 
from that of tensor product coefficients. Given such a choice the conjecture boils down 
to a simple relation for the minimum level of a coupling in terms of the corresponding 
minimum level of the elementary couplings appearing in its decomposition. Writing the 
decomposition in the form 

A <g> p, <g> v = Yl Ei* 

i 

the precise relation is 

k = ^biei (6.5) 

i 

where is the minimum level for the elementary coupling Ei. 

The generating function for su(S) fusion coefficients obtained from this method has 
been proven to be equivalent to the combinatorial algorithm of Cummins [22]. This con- 
stitutes a proof of the relation (6.5) for su(3). The correct choice of forbidden coupling 
turns out to be E1E3E5. 

Conversely, granting (6.5), it is simple to find which coupling must be forbidden. 
Indeed, since EiE 3 E 5 = E 7 E 8 = (1, 1)(1, 1)(1, 1), let us consider this specific fusion rule 
(c.f. example 8). We know that this coupling comes in two copies, with ko = 2,3. The 
corresponding BZ triangles together with their decomposition into elementary couplings 
are 

E2E4EQ E1E3E5 = E-jEs 

1 . . 

11 ( 6 - 6 ) 

11 
10 1 110 
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For su(3) ej = 1 for all i = 1, 2, . . . , 8, and therefore E1E3E5 must be forbidden in order 
to obtain the correct result ko = 2 for the second triangle. 

Having determined the forbidden coupling appropriate for the description of fusion 
rules, one can decompose uniquely any coupling into elementary couplings and read off ko 
from (6.5). By forbidding E1E3E5 it is straightforward to obtain the general and unique 
decomposition of a BZ triangle into a product of the Ei's. For instance, assuming that 
a 2 < min(a5, as), one finds the following decomposition for (6.2): 

E a 1E a B -a 2E a4 E a 3 -a 6E ar E a 2E a 6 ? (a 2 < min(a 5 , a 8 ) ) 

from which it follows that 

ko = a 4 + v\ + h> 2 , (a 2 < min(a 5 , a 8 )) (6.7) 

in agreement with (6.4) .i This proves that the strong depth rule with a c d is correct for 
su(3). 

Example 13: 

I 



The triangle E2E4E0 (k = 3) in (6.6) is associated with the standard tableau 

})) with a c e = 2 and thus ko = 
with a a = and thus ko = 2. 



2_ 

(c.f. the relation between (6.2) and (6.3)) with a c e = 2 and thus ko = 3, while the triangle 
E 7 E S (ko = 2) is associated with 



7. Conclusion. 

We have reformulated in a precise way the depth rule of Gepner and Witten [1], 
and shown how it can be used in principle to calculate fusion coefficients. However, the 

4 An equivalent result has been obtained in [23]. However, their approach is different since 
they started from the Kac- Walton formula, which relates fusion coefficients to tensor product co- 
efficients using the action of the affine Weyl group. Furthermore the expression for ko given by Lu 
is different from the one above. This is because his triangle is related to ours by a Schiitzenberger 
involution (t\t2t\) of the associated standard tableau. Explicitly, Lu's result is 

k (tit 2 ti~BZ) = max[min(a 3 , a 8 ) + /ii + ^2, min(a 2 , a e ) + v\ + ^2,min(a 5 , a 9 ) + Ai + A 2 ]. 



20 



cornerstone of the depth rule is the depth itself and we have argued that its calculation 
is problematic. More precisely, the depth rule is a constraint imposed on tensor product 
coefficients. But the basis convenient for calculating tensor product coefficients is not 
convenient for the calculation of the depth. 

To carry on, we compute the depth in a special (q = 0) limit of the quantum version of 
the finite part of the affine algebra under consideration, and freely extrapolate the result. 
Unfortunately, for su(N) (N > 3) it does not give the correct results. However we stress 
that counterexamples have been found only for moderately high levels and multiplicities. 
This indicates that the depth rule with crystal depth gives rather accurate lower boundsi 
for the minimum levels at which couplings first appear. On the other hand, for su(3) this 
procedure gives exact results. It (namely, (6.3) with ao given by (4.4)) turns out to be 
a rather efficient and simple algorithm for the calculation of fusion coefficients. A little 
puzzle remains, however: why does it work for sw(3), and why is it so accurate in general? 

As for sw(3), it seems that the success of the present approach is linked to the fact that 
all elementary couplings have level one, which implies that the decomposition of a coupling 
into elementary couplings can be done column by column in Young tableaux. In this vein, 
the failure of the depth rule with crystal depth charge for su(N > 3) could be explained 
by the existence of elementary couplings with minimum level higher than 1. However, 
thanks to syzygies, they can most of the time be eliminated from the decomposition of 
a coupling into elementary couplings. This should account for the reasonable success of 
the crystalline version of the depth rule. Finally, for sw(3), we noticed that BZ triangles 
provided a very convenient ground for the analysis of generating functions. We will report 
elsewhere on the extension of this approach to su(N > 3). 

We have argued that in the absence of an appropriate basis, generic fusion coeffi- 
cients cannot be easily computed from the depth rule (which is based on the consideration 
of three-point functions). But we should mention another way by which they could be 
extracted solely from the analysis of singular vectors. Indeed, among the Kac-Moody 
singular vectors are those that may be interpreted as mixed Virasoro-Kac-Moody singu- 
lar vectors, given the Sugawara construction. These mixed singular vectors lead to the 
Knizhnik-Zamolodchikov equation [10]. From this equation, one can evaluate the four- 
point functions and extract from these the three-point functions, which then yield the 

5 We have no proof that it is actually a lower bound, but the statement is motivated by the 
very large number of examples we have worked out. 
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fusion coefficients. However, to our knowledge, no such calculations have been performed 
for three-point functions such that N\n P > 1. 

On the other hand, the depth rule has been used in [6] to compute arbitrary fu- 
sion coefficients using a Schubert type calculus.i This does not contradict our conclusion 
since the depth rule is used there only to calculate fusion coefficients involving at least 
one fundamental representation (i.e. the truncated Pieri's formula). For a fundamental 
representation, the depth charge is easily evaluated and it agrees with its crystal version. 
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Appendix A. Tensor products with standard tableaux: words, bumping and 
the Robinson-Schensted correspondence. 

In this appendix we describe a more general way of representing states in a given 
coupling and of calculating tensor product coefficients. To a standard tableau one may 
associate a word, which is the sequence of numbers appearing in the tableau, read column 
after column from bottom to top and from left to right. 

Example 14: 

To the standard tableau 



1 


1 


2 


3 


5 


3 


4 


4 






5 











It should be stressed that this approach is conceptually very different from the method 
of generating functions [5]. Starting from the fact that any irreducible representation can be 
obtained from an appropriate multiple product of fundamental representations, it is shown in [6] 
that this description can be extended to the affine case provided we impose some level dependent 
constraints, called syzygies, on the ordinary product of the fundamental representations. On 
the other hand, for the method of generating functions, the basic variables are the elementary 
couplings and the syzygies there refer to relations between products of elementary couplings, 
already present for the calculation of tensor products. 
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we associate the word [531414235]. 

On the other hand, given a word one can reconstruct the corresponding standard 
tableau by the bumping method. One places the leftmost number in a box, the upper left 
box of the tableau being constructed. For the next number one proceeds as follows: If 
it is greater than or equal to the number already placed in the first box, it is put in the 
second box of the first row; otherwise it takes the place of the first number in the first 
box, the latter being bumped in the first box of the second row. Proceeding in this way 
for all numbers of the word and allowing bumping in every row, one easily reconstructs 
the corresponding standard tableau. This relation between standard tableaux and words 
is known as the Robinson- Schensted correspondence [26]. 

Example 15: 

Here is the step-by-step reconstruction of the tableau associated with the word [3121]: 



[3]^ [121] = 

(This example illustrates the fact that words are invariant with respect to some rear- 
rangements in the order of the entries: the words [3121] and [3211] are manifestly equiva- 
lent.) 

Any intermediate step in the reconstruction of a standard tableau from a word can 
be viewed as the insertion of a word into a standard tableau. The representation of a 
standard tableau by a word, followed by the insertion of this word into another standard 
tableau defines a 'fusion' procedure for two standard tableaux which exactly reproduces the 
Littlewood-Richardson rule for calculating tensor products coefficients. This construction, 
due to Thomas [27], provided the first combinatorial proof of the Littlewood-Richardson 
rule. More explicitly, in order to obtain all possible v which could appear in the product 
A ® ft one proceeds as follows. Consider all the states T(A^) (resp. T(/2^))@ in the A 
(resp. p) representation. By fusing all possible pairs of standard tableaux like T(A^) <— 
W(T(p, , (- ;■))), one fills up all the highest weight representations V which appear in the 
product A ® \x. Notice that the symmetry of the tensor product with respect to the 

7 In our notation T(A) is the standard tableau associated with the state |A), and W(T(X)) is 
the corresponding word. 



I 

3 



[21] 



[1] 



1 


1 


to 




3 
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interchange of A and n is reflected by the existence of a dual word W* and a right insertion 
operation such that [27] 



W*(T(\')) -> T{p!) = T(X') <- W(T(//')) 



Example 16: 

Fusion of the ST 

is done as follows: 



1 


2 


3 


2 


4 





and 



1 


4 


3 




4 





1 


2 


3 


2 


4 





[4314] = 



1 


1 


3 


3 


4 


2 


2 


4 





which shows that the representation (2, 2, 1) appears in the su(4) tensor product (1, 2, 0) ® 
(1,0,1). 

The generalization of the Thomas construction to the case of a triple product 
is straightforward: one simply has a chain of insertions of the type [T <— W(T')] <— W(T"). 

Example 17: 

Let us identify the two triplets of states contributing to the stt(3) tensor product 
(1, 1) <g> (1, 1) <g> (1, 1) D (0, 0). These are 



and 



Indeed, in the first case one has 



[312] = 



1 


2 


3 


2 


3 




3 







followed by 



1 


to 


3 


2 


3 




3 







[211] = 



1 


1 


1 


2 


2 


2 


3 


3 


3 



= (0,0) 



In this example, associativity is easily verified, i.e. 

[T <- W{T')] <- W(T") = T <- W^T' <- W(r /; )] 
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Example 18: 

With this method one can redo example 9. In order appropriate to left insertion, the 
contributing tableaux are 



2 


3 


3 


4 




1 


1 


1 


1 


3 


4 


4 




, x , 


2 


2 


2 




4 










3 









where X stands for five tableaux given in Ex. 9. Since a c e for the two other tableaux above 
is easily checked to be 4, one reobtains the set {4, 5, 5, 5, 6} for k . 
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